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We analyze the dispersion correction to elastic parity violating electron-proton scattering due to 
'yZ exchange. In particular, we explore the theoretical uncertainties associated with modeling contri- 
butions of hadronic intermediate states. Taking into account constraints from low- and high-energy, 
parity-conserving electroproduction measurements, choosing different models for contributions from 
the non-resonant processes, and performing the corresponding flavor rotations to obtain the elec- 
troweak amplitude, we arrive at an estimate of the uncertainty in the total contribution to the 
parity- violating asymmetry. At the kinematics of the Q-Weak experiment, we obtain a correction to 
the asymmetry equivalent to a shift in the proton weak charge of (0.0054 ± 0.0020). This should be 
compared to the value of the proton's weak charge of Cf^ = 0.0713 ± 0.0008 that includes Standard 
Model contributions at tree level and one-loop radiative corrections. Therefore, we obtain a new 
Standard Model prediction for the parity- violating asymmetry in the kinematics of the Q-Weak ex- 
periment of (0.0767 ± 0.0008 ± 0.0020 7 z). The latter error leads to a relative uncertainty of 2.8% in 
the determination of the proton's weak charge, and is dominated by the uncertainty in the isospin 
structure of the inclusive cross section. We argue that future parity-violating inelastic ep asym- 
metry measurements at low-to-moderate Q 2 and W 2 could be exploited to reduce the uncertainty 
associated with the dispersion correction. Because the corresponding shift and error bar decrease 
monotonically with decreasing beam energy, a determination of the proton's weak charge with a 
lower-energy experiment or measurements of "isotope ratios" in atomic parity- violation could pro- 
vide a useful cross check on any implications for physics beyond the Standard Model derived from 
the Q-Weak measurement. 

PACS numbers: 



INTRODUCTION 



Precise measurements of low-energy observables can 
provide powerful probes of physics beyond the Stan- 
dard Model that complement high energy collider stud- 
ies [2 [2] . In particular, measurements of parity- violating 
(PV) observables in atomic physics and electron scatter- 
ing have provided key tests of the neutral weak current 
sector of the Standard Model and constrained possible 
new physics in this sector [21 0HS]. In this work, we con- 
sider parity-violating (PV) elastic scattering of longitu- 
dinally polarized electrons from hydrogen, which is the 
subject of the Q-Weak experiment at the Jefferson Lab 
(JLab) [7 a . This experiment draws on a rich history of 
parity-violating electron scattering (PVES) at various fa- 
cilities and aims to provide the most precise determina- 
tion of Qyyj t ne weak charge of the proton, ever made. 

In PVES, the weak charge is operationally defined 
through the forward scattering limit of the PV asym- 
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metry: 
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where the ratio of response functions is defined below. 
Here and in the rest of the article, Gf denotes the Fermi 
constant, as taken from the muon lifetime (often denoted 
by G^). The weak charge - defined as a static property 
of the proton - is then the leading term the expansion of 
the ratio W PV /W EM in powers of t = -q 2 : 
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where the reason for specifying zero beam energy E 
will become apparent below. In the one-boson exchange 
(OBE) approximation, the weak charge is just given by 
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where the C\ q characterize the effective four-fermion 
parity-violating electron-quark interaction 
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In the Standard Model, it is possible to make pre- 
cise predictions for the Ci q , including the effects of 0(a) 
electroweak radiative corrections [31 El El - These correc- 
tions include the effects of one-loop contributions to the 
gauge boson and fermion propagators and gauge boson- 
fermion vertices. Ultraviolet (UV)divergences are re- 
moved through renormalization, and in what follows we 
will use the modified minimal subtraction (MS ) scheme 
for doing so. 

Additional, UV-finite corrections arise from the two- 
boson exchanges ("box graphs"): ZZ, W + W~, 7Z and 
77. Those involving two heavy vector bosons are dom- 
inated by loop momenta of order Mz and are properly 
included in the radiatively-corrected C\ q coefficients. On 
the other hand, the box graph corrections involving one 
or more photons are sensitive to low-momentum scales 
where target-dependent hadronic structure effects may 
be significant. In what follows, we focus on the Z7 box 
correction. For a review of recent work on the 77 correc- 
tions, see Ref. [TO] , 

Recently, the 7Z box graph contribution has been the 
subject of renewed scrutiny. In Refs. [SHHE], the short- 
distance part of this correction was computed, confirming 
the earlier computation of Ref. [5]. It carries a logarith- 
mic dependence on the hadronic scale, Ahad, with the 
latter requiring the presence of a "low energy constant" 
C7z(Ahad) to yield a result independent of the hadronic 
matching scale. The authors of Ref. [31 [5] assigned a 
generous error to C 7 z(Ah a d) associated with the difficult- 
to-computc long-distance hadronic effects. 

The authors of Ref. [TT] subsequently observed that 
there exists an additional contribution from the 7.Z box 
graph that grows with the electron beam energy and that 
is independent of the hadronic cutoff parameter 1 . Given 
the energy-dependence of this "dispersion correction" , it 
is more appropriate to consider it as a new term in the 
PV asymmetry than as a contribution to the weak charge 
that is nominally a static property of the proton. Nev- 
ertheless, in the forward limit of Eq. ([2]), its effect is to 
shift the apparent value of Q v w - Moreover, unlike the 
short-distance and C 7 z(Ah a d) terms that are suppressed 
by 1 — 4 sin 2 9w ~ 0.07, the energy-dependent correc- 
tion is not accidentally suppressed. For the energy of the 
Q-Weak experiment, the authors of Ref. [TT] estimated 
that the correction was several percent, raising the pos- 
sibility that the estimated theoretical uncertainty in the 
PV asymmetry could be larger than given in Refs. [3j[6]. 

A follow-up study [H] repeated the computation of 
Refs. [TT][TB] using a somewhat different hadronic model 
framework and drawing upon recent structure function 
measurements carried out at the Jefferson Laboratory. 
These authors argued that the expressions used in Ref. 



For related work considering the effects of the ■yZ box graph away 
from the forward limit — relevant to the strange quark form factor 
determinations - see Refs. I12H14| . 



[TT] contained numerical errors but nonetheless obtained 
a quantitatively similar result for the size of the correc- 
tion. An estimate of the uncertainty in the correction 
was also provided, suggesting that the theoretical uncer- 
tainty associated with the energy-dependent term is well 
below the uncertainty quoted in Refs. [34 E]. Recently, 
another study of this correction was reported in Ref. [17] . 
The latter work employed yet another parametrization of 
virtual photoabsorption data from Jefferson Lab, and a 
different treatment of the isospin structure and of the 
uncertainty was applied. The results is consistent with 
that of Ref. [TS] with an error bar that is also smaller than 
that of Refs. [3J E] . We will review these works in greater 
detail below. For the moment, we display in Table [TJ the 
results of the previously mentioned studies along with 
the results of this work. While all of the recent results 
(ours and Refs. [13 H2]) are consistent within quoted er- 
ror bars, we obtain a larger uncertainty by roughly a fac- 
tor of two. As we discuss below, this larger theory uncer- 
tainty results from taking into account hadronic model- 
dependence in computing the jZ dispersion correction. 



Ref. [TT] 


Ref. [15] 


Ref. [17] 


This work 


(3±3)10 -3 


(4.7tJ;I)10- 3 


(5.7±0.9)10 -3 


(5.4 ± 2.0)10" 3 



TABLE I: Estimates for the dispersion correction Re D 7 ^ ob- 
tained in various works, as indicated in the Table. Originally, 
Gorchtein and Horowitz in [TT] only quoted the value of Re 
□ 7 # ~ 6%, as calculated relative to 1 — 4 sin 2 8w(0) « 0.05. 
This corresponds to the number given in Table. 



Obtaining a robust theoretical prediction for A in 
the Standard Model is essential for the proper interpreta- 
tion of the asymmetry in terms of possible contributions 
from physics beyond the Standard Model. In light of the 
recent history and disagreements in the literature on the 
question of the 7Z box correction, we revisit here the 
computations of Refs. [TTJ [T5Hl7| . Our goal is three- 
fold. First, we seek to clarify the apparent disagreements 
about the numerical factors in the analytic expressions 
for the energy-dependent part of the "fZ correction. Sec- 
ond, we attempt to provide an estimate of the theoretical 
uncertainty associated with hadronic modeling required 
for its computation. While the study of Ref. [15] in- 
cluded an uncertainty associated with the experimental 
data used as input for the calculation, no estimate of the 
theoretical error related to the choice of model framework 
was given. Finally, we discuss additional experimental in- 
put that would be useful to improve the reliability of the 
calculated correction. 

The remainder of our treatment of these points is orga- 
nized as follows. Section II outlines the elastic electron- 
nucleon scattering kinematics and observables that are 
analyzed to one-loop order. In Section III, we derive a 
forward dispersion relation for the dispersion corrections. 
In Section IV, we discuss the input in these sum rules, 
perform an isospin decomposition of the inclusive electro- 
production data and isospin-rotate these data in order to 
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obtain the inclusive parity violating data. We combine 
different data sets to obtain an estimate of the uncer- 
tainty associated with such rotation in the flavor space. 
Detailed discussion of the isospin rotation of the resonant 
contributions is reported in Appendix A. In Section V, 
we present our results for the dispersion correction D 7 z 
and the respective theory uncertainty at the kinematics 
of the QWEAK experiment. Section VI is dedicated to 
the study of the t-dependence of the dispersion correction 
that is important for translating the value obtained from 
dispersion relation in the exact forward direction to the 
experimental kinematics. In Section VII, we compare the 
existing calculations of the energy-dependent dispersion 
correction to the weak charge of the proton in detail. 



We close the article with a short summary in Section 
VIII. 



II. PVES IN THE FORWARD SCATTERING 
REGIME 

We consider elastic scattering of massless electrons off 
a nucleon, e(k) + N(p) — > e(k') + N(p'), in presence of 
parity violation (and in absence of CP- violation). The 
scattering amplitude T can be cast in the following form 
involving six scalar amplitudes fi(y, t), i — 1, 2, 6, 
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where only electromagnetic and weak neutral currents are 
considered. Gf stands for the Fermi constant, as taken 
from the muon lifetime, according to the MS scheme. The 
amplitudes fi.2,3 are parity conserving (PC), and /4,5,6 
are explicitly parity violating (PV). Above, k(k') stands 
for the initial (final) electron momenta, and pip 1 ) for the 
initial (final) nucleon momenta, respectively, and M de- 
notes the mass of the nucleon (we take M n « M p = M) . 
All six amplitudes are functions of energy v = (with 
K = an( j p — 2±£-) and the elastic momentum 

transfer is t = A 2 < 0, with A — k — k' = p' — p. At tree 
level (one boson exchange, OBE) and to leading order 
in Gf and a em , the amplitudes fi reduce to the elec- 
tromagnetic and weak form factors of the nucleon (the 
index N takes values p, n denoting proton and neutron, 
respectively) , 
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with t = 



In absence of radiative corrections, these 
amplitudes reduce to the electroweak Sachs form factors 
G'^ Z j. In terms of these generalized form factors, the 
unpolarized cross section on a nucleon target N can be 
written as 
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with 9 the electron Lab scattering angle, -E(-E') the in- 
coming (outgoing) electron Lab energy, and e = (1 + 
2(1 +r) tan 2 |) _1 the virtual photon longitudinal polar- 
The reduced cross section a^, up to 



ization parameter 
and including terms of order a 



reads 
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Above, gy = —(1 — 4sin 2 9w) and g e A — 1. Radiative 
corrections induce terms 6 fa ~ a em , leading generically 
to fi = f? BE it) + Sfiiu,t). We denote the usual Dirac 
(Pauli) form factors by J~i2j respectively, and the nu- 
cleon axial form factor at tree level by G\ N . Similarly, 
•F\2 stand for the form factors describing the vector cou- 



pling of the Z to the nucleon. 
ventional combinations, 

f N + 



One introduces the con- 
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fi 

fx - rf 2 
1 

a, 



In what follows, we will concentrate on the case of 
electron-proton scattering. Therefore, we will under- 
stand N = p everywhere and suppress the index N in 
all expressions, unless explicitly stated otherwise. 

The parity violating asymmetry is defined in Eq. ([T]) 
with the ratio of the response functions is given by 
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Here, <r± are the cross sections for positive and negative 
helicity electrons, and e' = yr(i + t)(1 — e 2 ). 

Since we are interested in very forward scattering an- 
gles sa 8° corresponding to the Q-Weak kinematics [7], 
thus r < 10 -3 , the expressions for the cross section and 
PV asymmetry can be fruther simplified. 

For the reduced cross section the leading contribution 
in Eq. ^ comes from the G\ term, and we obtain 



1 



OR 



(77)^(1 + T8% in + 2Re^ c ) + 2ReO rr . (11) 



The three distinct corrections quoted above are defined 
as follows: <5£ in is a kinematic correction that arises at 



tree level due to the magnetic part and other sublead- 
ing kinematic effects of order r, (1 — e), that do not con- 
tain 0(a em ) effects; 5 RC stands for order 0(a em ) cor- 
rections that are energy-independent (such as vacuum 
polarization, self energy and vertex corrections) ; finally, 
□ 77 denotes the two-photon exchange correction that is 
an energy-dependent 0(a em ) correction. 

Similarly, for the PV asymmetry the leading order con- 
tribution in Eq. (10 1 originates from the G§*Gg term. 



As discussed in Ref. [6J, the Standard Model predic- 
tion for the PV asymmetry in the forward regime can be 
expressed as 
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where 9w (0) is the running weak mixing angle in the MS 
scheme at zero momentum transfer [3]. The correction Ap 
is a universal radiative correction to the relative normal- 
ization of the neutral and charged current amplitudes; 
the A e and A' e give, respectively, non-universal correc- 
tions to the axial vector Zee and jee couplings; the D vv 
for V — W, Z, 7 give the non-universal box graph correc- 
tions; and the "+ • • • " indicate terms that vanish with 
higher powers of t in the forward limit, such as those 
arising from the magnetic and strange quark form fac- 
tors and the two-photon dispersion correction, D 77 . The 
weak charge of the proton, considered as a static prop- 
erty, is given by the quantity in the squark brackets in 
the zero-energy limit. 

Within the radiative corrections, the TBE effects are 
separated explicitly. This is done because the TBE cor- 
rections, unlike other corrections in the above equation, 
are in general v and i-dependent. In particular, the v (or 
e) dependence of the 77-box is believed to be responsible 
for the discrepancy between the Roscnbluth and polar- 
ization transfer data for G~^,/G M [15] . It should be noted 
that in the exact forward direction D 77 vanishes as a con- 
sequence of electromagnetic gauge invariance. 

The WW and ZZ-box diagrams were first considered 
in [3], and subsequently investigated in Refs. [5], [T5] . The 
contribution from □ww in particular is relatively large. 
Both corrections are ^-independent at any hadronic en- 
ergy scale since they are dominated by exchange of hard 
momenta in the loop ~ Myy, Mz- Higher-order perturba- 



tive QCD corrections to Dvkw an d Ozz were computed 
in Ref. [6], and the overall theoretical uncertainty associ- 
ated with these contributions is well below the expected 
uncertainty of the QWEAK experiment. 

In contrast to □ww and Dzz, 7 z receives substantial 
contributions from loop momenta at all scales. For the 
electron energy-independent contribution, this situation 
leads to the presence of a large logarithm lnM^/Ahad 
where A had is a typical hadronic scaled [HI HHj- Since 
the asymmetry must be independent of the latter, O-yz 
includes also a "low-energy constant" C 7 z(Ahad) whose 
hadronic scale dependence compensates for that appear- 
ing in the logarithm. An analogous 14^7 box correction 
enters the vector current contribution to neutron and nu- 
clear /3-decay. Importantly for the PV asymmetry, these 
energy-independent 7Z box contributions are suppressed 
by 1 — 4 sin 2 9\y, thereby suppressing the associated the- 
oretical uncertainty. 

In Ref. [llj , the 7^-box contribution was re-examined 
in the framework of dispersion relations and it was found 
that it possesses a considerable energy dependence, so 
that at energies in the GeV range its value can differ 
significantly from that found at zero energy. Moreover, 
the energy-dependent correction contains a term that is 
not 1 — 4sin 2 #^ suppressed, so the theoretical uncer- 
tainty associated with hadronic-scale contributions is po- 
tentially more significant. This energy dependence comes 
through contributions from hadronic energy range inside 
the loop that cannot be calculated reliably using pertur- 
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bative techniques. 

At present, a complete first principles computation is 
not feasible, forcing one to rely on hadronic modeling. 
For a proper interpretation of the PV asymmetry, it is 
thus important to investigate the theoretical hadronic 
model uncertainty. The remainder of the paper is de- 
voted to this task. In so doing, we will attempt to re- 
duce this model uncertainty by relating - wherever pos- 
sible - contributions from hadronic intermediate states 
to experimental parity-conserving electroproduction data 
through the use of a dispersion relation and isospin ro- 
tation. As a corollary, we will also identify future ex- 
perimental measurements, such as those of the parity- 
violating inelastic asymmetry in the regime of moderate 
Q 2 and W, that could be helpful in reducing the theo- 
retical uncertainty. 



III. DISPERSION CORRECTIONS 

To calculate the real part of the 7Z direct and crossed 
box diagrams showed in Fig. [T] we follow [11] and adopt 
a dispersion relation formalism. We start with the calcu- 
lation of the imaginary part of the direct box (the crossed 
box contribution to the real part will be calculated using 
crossing) , 



ImT, 



7 z 



d 3 h 



b fJ,U 



V2 (2^) 3 J 2E l Q 2 (l + Q 2 /M 2 )' 
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where Q 2 — — {k — fci) 2 denotes the virtuality of the 
exchanged photon and Z (in the forward direction they 
carry exactly the same Q 2 ), and we explicitly set the 
intermediate electron on-shell. In the center of mass of 
the (initial) electron and proton, one has E\ — , 
with s the full cm. energy squared and W the invariant 
mass of the intermediate hadronic state. Note that for 
on-shell intermediate states, the exchanged bosons are 
always spacelike. 





FIG. 1: Direct and crossed diagrams for 7Z-exchcange. 
Dashed lines correspond to an exchange of a Z-boson, and 
wavy lines to an exchange of a photon. The blob stands for 
an inclusive sum over intermediate hadronic states. 



The leptonic tensor is given by 



u(k')j l/ jt 1 j t ,(g'y + g e A j 5 )u(k). 



(14) 



We next turn to the lower part of the diagrams in Fig. 
[T] The blobs stand for an inclusive sum over all possible 
hadronic intermediate states, starting from the ground 
state (i.e., the nucleon itself) and on to a sum over the 
whole nucleon photoabsorption spectrum. The case of 
the elastic hadronic intermediate state was considered in 
20 . Here, we concentrate on the inelastic contribution. 
Such contributions arise from the absorption of a photon 
(weak boson). In electrodynamics, for a given material, 
the relation between its refraction coefficient and the de- 
pendence of the latter on the photon frequency (i.e., dis- 
persion) on one hand, and the photoabsorption spectrum 
of that material on the other hand, is historically called 
a dispersion relation. It is exactly this dependence of 
the forward scattering amplitude fi(i/,0) (see Eq. (5J)) 
on the energy that arises from its relation to the elec- 
troweak 7(Z)-absorption spectrum that is the scope of 
an investigation in this work. This explains the origin 
of the term "dispersion correction" used for the inelastic 
contributions to the ^yZ-boK correction. 

In the forward direction, the imaginary part of the dou- 
bly virtual "Compton scattering" (j*p — > Z*p) ampli- 
tude is given in terms of the interference structure func- 
tions 3(^5 Q 2 ), with x = the Bjorken variable. 
Making use of gauge invariance of the leptonic tensor, 
we have 



F, 



-yZ 



ic ^a0 Palp pjZ 



Pq 



(15) 



Contracting the two tensors, one obtains after a little 
algebra two contributions that are due respectively to the 
axial and vector couplings of the Z to the electron, 
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where the imaginary parts ImD will appear in a disper- 



sion relation for the real parts in Eq. (20) below. The 
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full correction is the sum of the two, 



ImD 7 z(z/) = lmD yZA (u) +ImD jZv (v). (17) 



In Eqs. (16 1, W 2 = (M + m^) 2 stands for the pion pro- 
duction threshold, and the Q 2 -integration is constrained 
below a maximum value 



Q 



2 

max 



(s- M 2 )(s~W 2 ) 



(18) 



as a condition of on-shell intermediate states for an imag- 
inary part calculation. Eq. ( 16 ) is in agreement with 



Refs. P21 HZ] . In particular, we confirm the correctness 
of the claim made in Ref. [TS] that in Ref. [TT] a factor 
of 2 was missing. 

In order to write down the dispersion relation for the 
function IZ] 7 z(i/), one should consider its behavior under 
crossing. We distinguish two contributions, O jZv and 
□ 7 z A that have different crossing behavior [TT] : 



□ 7 z A (-i>) = -n jZA (v) 
u lZv {-v) = +n lZv {v) 



(19) 



Correspondingly, the two contributions obey dispersion 
relations of two different forms, 



ReD 7 z A {») 
ReD 7Zv (v) 



2v 
2 

TT 



v' 2 — v 
v'dv 
v 12 — V 



2 ImD 7ZA (z/) 
2 lmU lZv { V ') (20) 



where the presence or absence of the factor of v' in the 
integrands follows from the behavior of the ImD under 
crossing symmetry. 



The result in Eq. ( 20 ) gives a model-independent rela- 



tion between the dispersion correction to the weak charge 
of the proton and the parity violating structure functions 
appearing in Eq. ( 16 1. This relation does not rely on any 



assumption, other than the neglect of higher order radia- 
tive corrections and the number of subtractions needed 
for convergence of the dispersion relation. The advan- 
tage for this formulation is that the F^ z are in principle 
measurable. However, in absence of any detailed parity 
violating inclusive electron scattering data, the input in 
the dispersion integral will depend on a model. In the fol- 
lowing, we will investigate the extent to which this model 
dependence can be constrained by existing or future ex- 
perimental data. 



IV. INPUT TO THE DISPERSION INTEGRAL 

In the previous section, the contribution of the forward 
hadronic tensor to the box diagram was considered. In 
this section, we will address the possibility of relating the 
interference hadronic tensor of Eq. (15) 



d 4 ze^(iV|T[jr m (z)J^(0)] \N) (21) 
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FIG. 2: The plane Q vs. W and kinematic regions corre- 
sponding roughly to various physical contributions. 



to the pure electromagnetic one, 

W% = \J d*ze*'(N\T [J v em {z)J» m m | AO- (22) 

Using unitarity, we rewrite these matrix elements as 
an inclusive sum over intermediate hadronic states, 



4 ze iqz 



"£(N\JZ m (z)\X)(X\J» Cv (0)\N) 



x 



x 



(N\r NC (z)\X}(X\JU0)\N) 



(23) 



and 



\ f d 4 ze^Y,( N \ J e m (z)\X)(X\J? m (0)\N) 
J x 

(24) 



respectively. We now proceed to investigate the possible 
relationships between the products of transition matrix 
elements appearing in each inclusive sum (23) and (24 1. 



Theoretically, calculating the full set of contributions 
to the inclusive sum represents a fundamental difficulty 
since in QCD, the basis for intermediate states X is infi- 
nite, and the matrix elements are non-perturbative. Un- 
der certain kinematic conditions, one can organize this 
basis into leading and subleading (kinematically sup- 
pressed) sub-sets. We depict this situation schematically 
in Fig. [2j where we show in the Q 2 — W 2 plane the 
approximate kinematic areas where various mechanisms 
dominate. At high energy and Q 2 , and finite Bjorken 
x, the leading set of states is X = q + X' (q denotes a 
quark), where to leading order in l/Q, X' is a spectator. 
Thus, in this regime the electromagnetic (weak) current 
directly probes a single quark within the nucleon, and 
gives access to the parton distribution functions (deep 
inelastic scattering, DIS in Fig. [2]). At high energy and 
Q 2 , and small x, however, the picture changes, as the 
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leading set is X = qq + N. In this regime, the pho- 
ton polarizes the QCD vacuum at the periphery of the 
hadron, and the resulting gg-pair forms a color dipole 
that interacts with the nucleon (diffractive DIS in Fig. 
[2]). This picture was first realized in the Vector Meson 
Dominance model (VDM) that capitalized on the fact 
that since vector mesons and the photon have the same 
quantum numbers, the latter can fluctuate into former 
[5TJ [52] . This simple model works quite well at low Q 2 
(VDM area in Fig. |2j. Such "hadron- like" behavior of 
a photon in scattering processes also results in the e.-m. 
data following the Regge behavior, as a function of W (re- 
spective Regge area in Fig. [2]). At higher values of Q 2 , 
rescattering effects in vector meson-nucleon scattering 
become increasingly important but can still be accounted 
for in what is called the "generalized VDM" (GVDM re- 
gion in Fig. [2]). At low energies, the relevant degrees of 
freedom are hadronic (that is, highly non-perturbative), 
X = N, irN, ttttN, N* , A etc. In this regime, the inelastic 
cross section is tyically dominated by resonances on top 
of a non- resonant background (Resonance area in Fig. [2]). 
The boundaries of each kinematic region are, of course, 
approximate. Their meaning is that the farther one de- 
parts from a kinematical region, to the lesser extent the 
respective mechanism works. Consequently, a large area 
on the W 2 — Q 2 plane, that overlaps with all the de- 
picted regions but not covering them completely is the 
so-called shadow region where none of the mechanisms 
can be considered as fully dominant. 

If data for the interference cross section existed 
throughout all these distinct regimes, we would not need 
to know details of any of the aforementioned models. In 
principle, such data could be obtained with measure- 
ments of the PV inelastic asymmetries in the various 
kinematic regimes shown in Fig. [2] At present, however, 
either no or very poor data on PV inelastic scattering 
exist. Consequently, we will instead pursue an alternate 
strategy, endeavoring to make use of extensive data sets 
for real and virtual photoabsorption that exist through 
vast kinematic region in energy and Q 2 . To that end, we 
will rely on models that adequately describe the photoab- 
sorption cross section in different regimes and for each 
attempt to establish relationships between the matrix el- 
ements (X\ Jg m \N) and (X\ J^ C \N) for each intermediate 
hadronic state \X) of definite isospin. We will approach 
this problem by extracting the electromagnetic matrix el- 
ements from inclusive e.-m. data, and then isospin-rotate 
every such matrix element. We begin with a brief review 
of the experimental situation and discuss various model 
descriptions. 



A. Real and virtual photoabsorption data 

We find that the dispersion integral for Ke\3 1 z is dom- 
inated by moderate values of W < 5 GeV and Q 2 < 3 
GeV 2 (see Fig. 
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in Section V). Consequently, we need 
to analyze in detail contributions from the resonance 



regime and portions of what we have called the VDM, 
GVDM, and regge regimes. Our goal will be to draw 
upon existing experimental data for inclusive and semi- 
inclusive electromagnetic data to infer the interfer- 
ence structure functions that appear in the dispersion 
integrals. To that end, we first summarize the experi- 
mental situation. 

• Real photoabsorption cross sections have been mea- 
sured from the pion threshold to very high energies 



• Virtual photoabsorption data: high precision data 
from the JLab E94-110 [29] and the preliminary 
data from the E00-002 [30] experiments are avail- 
able in the resonance region; in the DIS region, 
we quote the data for the DIS structure function 
F 2 from SLAC NMC Collaboration [31], FNAL 
E665 collaborations [32] and DESY HI Collabo- 
ration [55] , 

While it is equally possible to use structure functions to 
describe resonance data, in the following we opt to use 
total photoabsorption cross sections with transverse or 
longitudinal (for virtual photons only) photon polariza- 
tion. These cross sections are unambiguously related to 
the electromagnetic structure functions, 



<tI p (W 2 ,Q 2 ) 



1 



2x W 2 - M 2 



W 2 - M 2 ' 
8ir 2 a 

" W 2 _ M 2 

2M 2 



(25) 



F^(^Q 2 )-Fr(x,Q 2 ) 



01 



with the usual Bjorken scaling variable x = w i_ M 2 + q2 ■ 
This choice is convenient because in what follows, we will 
address transitions between helicity states of the nucleon 
and resonances, and it is preferrable to work with matrix 
elements of the electromagnetic current with definite hc- 
licities. As is evident from Eq. ( p5| , the two helicity 
states are mixed in F 2 . Similar relations hold between 
the interference cross sections aj.£ p (W 2 , Q 2 ) and inter- 
ference structure functions F± 2 (x, Q 2 ). Note that the 
definition of the transverse and longitudinal polarizations 
of the photon and the Z-boson are identical since in both 
cases they are fixed by the lepton kinematics of the reac- 
tion e +p ->• e' + X. 

Real photoabsorption data exhibit the following gen- 
eral features: i) a resonance structure on top of ii) a 
smooth non-resonant background between the thresh- 
old of pion production and W ~ 2 — 2.5 GeV, and iii) 
Regge behavior at high values of W with the cross sec- 



tion that grows slowly with energy, a^. 



{W 



2\a r 



with tip ~ 1.095 the parameter of the pomeron. 



?Z(w 2 ) - URcs 

_7P 



a 



{W 2 ) + 
(W 2 ) 



- {w 2 ) a -\ 



(26) 



where a = ap,otf 2 , etc. stand for pomeron and Regge 
trajectories. In this work, the most recent fit in terms of 
two trajectories (pomeron plus / 2 ) is used [32] 



Tikgd 



= fi 



thr 



W 2 

(145.0 ± 2.0) m&I^ 



-0.5 



,'W 2 

(63.5 ± 0.9) fib ( ^ 



(27) 



with parameter of the pomeron e = 0.097 ± 0.002. The 
threshold factor / t hr is necessary to make the continua- 
tion of the Regge fit into the resonance region meaningful. 
In this work, we take it in the same form as in 1351 



/thr = 1 - exp 



W 2 - (M + m») 5 
M 2 



(28) 



For virtual photons in the range of W 2 , Q 2 of inter- 
est here, the picture remains the same, with the Q 2 - 
dependence of the resonance contributions described by 
the form factors measured for a number of resonances, at 
least in certain channels. We will next specify two mod- 
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FIG. 3: (Color online) World data on total photoabsorption 
[23H28] (see [22] for the complete list) compared to the two 
models described in the text. The experimental errors are not 
shown. 

els that provide a smooth extrapolation between the real 
photoabsorption data and the virtual photoabsorption 
data and that can to certain extent be used to describe 
data all the way up into the diffractive DIS region. The 
two models differ in the form of the Q 2 -dependence of 
the background contribution: 

• Model I: The model used in pj] utilized the res- 
onance parameters obtained in [35 and the non- 
resonant Regge contribution from [33] that was 
fitted to the real photoabsorption data at high 
energies. The Q 2 -dependence of the high-energy 
part was taken from the hybrid GVD/color dipole 
(CDP) approach of Ref. [37]. For the estimates of 



Parameter 


Ref. [3S] 


Model I 


Model II 


Cl 


2.124 


2.24 


2.2 


C2 


2.569 


2.73 


2.73 


C 4 


0.064 


0.155 


0.155 


C5 


0.549 


0.549 


0.7 


C6 


1.914 


1.914 


2.5 


C7 


1.0 


1.0 


1.5 


A 7 T (0) 


3.419 


5. 


5. 


Al(0) 


11. 


15. 


15. 



TABLE II: The list of the resonance parameters and their 
values for Model I and Model II, as compared to the original 
fit of Ref. [38] (see Table III of that reference) . The notation 
of 38 was kept. 



[TTj . a simple dipole model with the dipole mass 
A w 1 GeV for all the transition resonance form 
factors was employed. Because it was found that 
this simple dipole form fails dramatically through- 
out the resonance region, we adopt the resonance 
part from [38 with a few parameters minimally ad- 
justed in order to fit the data with the background 
of a different form, rather the one used in [35] orig- 
inally. We list those parameters and the respective 
changes in Table [TT] 

• Model II: To test the sensitivity of our calcula- 
tions to the specific model, we use another form 
of the background from the "naive" GVD model 
of Ref. [35] (cf. Eqs. (3,4) of that Ref.), and 
we add the resonance contributions from [38] on 
top of that. Again, some resonance parameters are 
slightly adjusted to the background, and all changes 
are quoted in Table [TT| 

In Fig. [3] we confront the two models with the total 
photoabsorption cross section. The Model I is shown by 
solid red lines, Model II by the dashed blue line. 

Figs. [4|6] display the comparison of the two models 
with the data for the differential cross section for inclusive 
electroproduction in the resonance region. Both models 
in general provide a good description of the data in the 
resonance region. The areas between the lower and up- 
per thin curves in each plot correspond to the range of 
values of the helicity amplitudes for the photoexcitation 
of each resonance included in Models I and II, as given 
by the PDG [3D]. It can be seen that the experimental 
data are always contained within these areas for W 2 < 4 
GeV 2 , even without including the experimental errors. 
At the same time, we note that just above the resonance 
region, in the limited range 4 GeV 2 < W 2 < 6 GeV 2 , 
and at moderate values of Q 2 , the background system- 
atically lacks strength. However, we stress that this lack 
of strength is observed only in very limited range of en- 
ergies, and the deficit is less than 20% which makes the 
impact of this effect on the dispersion correction small. 
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FIG. 4: (Color online) Differential cross section data in the 
resonance region from [30] are shown in comparison with the 
two models. The experimental errors are not shown. Thick 
solid line is the result of Model I, and thick dashed line is the 
result of Model II. Thin solid lines show the error bar due to 
the uncertanties in helicity amplitudes for the photoexcitation 
of the resonances on the proton, according to [40] . 



We next turn to the deep inelastic (DIS) data. For DIS, 
a natural choice would be to use the PDF parametriza- 
tions from MRST or CTEQ, DGLAP-evolved to the nec- 
essary value of Q 2 . However, this is only applicable at 
large enough Q 2 , and extrapolating them below Q 2 = 1 
GeV 2 introduces additional systematic error. In Figs. 



FIG. 5: (Color online) Differential cross section data in the 
resonance region from [30j are shown in comparison with 
Model I and II. Notation as in Fig. [4] 
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11 



the naive GVD model of Ref. [35] (Model II) is 
shown along with the GVD/CDP model of [37J (Model I). 
One can see that while the GVD / CDP model reproduces 
the data in a wide range of x, Q 2 , the naive GVD model 
overshoots the data at large x starting at moderate Q 2 , 
and underestimates the low-x behavior for all Q 2 . One 
needs to keep in mind, however, that both models work 
reasonably well at moderate Q 2 and large x which give 
the main contributions to the dispersion correction. 
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FIG. 6: (Color online) Differential cross section data in the 
resonance region from f2U] are shown in comparison with 
Model I and II. Notation as in Fig. [4] 
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FIG. 7: (Color online) Differential cross section data in the 
resonance region from [29] are shown in comparison with 
Model I and II. Notation as in Fig. [4] 



The following comment is in order here. The authors 
of Ref. [TS] argued that our description of the data is 
unsatisfactory not only in the resonance region but also 
beyond (cf. Fig. 1 of Qjjj). While the model of the 
resonance form factors of Ref. [IT] was definitely not ac- 
curate (one of the instances on which we improve that 
calculation in the present work), the model for the back- 
ground in [TTJ is exactly the same as that of Model I here. 
We believe that Figs. |4fTT1 presented in this section pro- 
vide abundant evidence of a satisfactory description of 
the experimental data by our phcnomcnological model. 
In view of this, we find it puzzling that Ref. [15] quotes 
a discrepancy of 40-50% at Q 2 as low as 0.6 GeV 2 just 
above the resonance region (cf. the upper left panel of 
Fig. 1 of that reference). 



B. Isospin rotation of the resonance contributions 

In Standard Model, the Z and 7 hadronic currents are 
related by means of a simple isospin rotation, 

JL = 9 /=0 ^=o + ? /=1 ^=i + « 

J&Cy = 9v = °Ji=o + 9v =1 JLi + 9 s yJ^ (29) 

with 

J% = S7**s . (30) 
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FIG. 8: (Color online) Differential cross section data in the 
resonance region from are shown in comparison with 
Model I and II. Notation as in Fig. [4] 
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The e.m. charges given by 

1 i=i 1 



1=0 



3^/2' q y/2' 
whereas the weak charges are 

1 4 



sir 



sfr 1 



ffv 



713 s Uw 



1 

V2 



(2 - As 2 e w ) 



-1 + ~s J 



(31) 



with s 2 9w being a shorthand for sin 2 By/ (for purposes 
of this argument). This isospin decomposition is used 



FIG. 9: (Color online) Differential cross section data in the 
resonance region from [25] are shown in comparison with 
Model I and II. Notation as in Fig. [4] 



to relate weak proton form factors to the proton and 
neutron electromagnetic form factors, 



W&c.vlP) = (1 - 4s 2 M(H J£n\p) - (n\JL\n) (32) 



where we neglected strangeness contributions that are 
generally small [4"T] . 
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FIG. 10: Comparison of the x-dependence of the DIS struc- 
ture function F2(x,Q 2 ) at fixed Q 2 and as a function of x, 
in GVD/CDP model of [37] (solid lines) and the naive GVD 
model of [39] (dashed lines) to the low- a; DIS data of HI Col- 
laboration [31]. The experimental errors are not shown. 
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TABLE III: Ratios yn with respective uncertainties for seven 
resonances. 



The above relation is valid for transitions to / = | 
resonances, as well: 



(X\J*\p) = (1 - 4s 2 9 w )(X\J? m \p) - (X|J£» (33) 



It is then straightforward to relate the contribution of 
a resonance R with isospin 1/2 to the interference -fZ 
cross section entering Eq. ( 16 ) to its contribution to the 



electromagnetic cross section: 



(p\ja m \R}(R\j»\p) 



(l-4s*6 w )\{R\J? m \p}\ 2 
(p\J? m \R)(R\J? m \n) 



(34) 



Consequently, for each resonance, we define two ratios 
describing the relative strength of its contribution to the 
7Z-interferencc cross sections erj^ R with respect to the 



purely electromagnetic ones cr^ p , L -. R as 



Cl /7 (Q 2 ) 



7 T,R 
U T,R 

'L,R 



L.R 



(35) 



In the Appendix A we discuss in detail the (in- 
dependence of these ratios, as well as the ratios of the 
longitudinal cross sections Cz/y- Basing on the discus- 
sion in Appendix A, we will use the value 

tz h {Q 2 ) = [l-4.s 2 M0)] -VR = const., (36) 

to rescale the contribution of a resonance R to both trans- 
verse and longitudinal cross section. Possible discrepan- 
cies (which, if known, are model-dependent) from this 
rule are accounted for by assigning a conservative uncer- 
tainty to the ratios £^ , . This is done by using the PDG 
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FIG. 11: Comparison of the Q 2 -dependence of the DIS struc- 
ture function F2(x,Q 2 ) at fixed x and as a function of Q 2 , 
in GVD/CDP model of [37J (solid lines) and the naive GVD 
model of [39] (dashed lines) to the DIS data of NMC Collab- 
oration [3T] and E665 Collaboration [21] where the :r-birming 
corresponds to that of NMC. The experimental errors are not 
shown. 



values and respective errors for the transition helicity am- 
plitudes. These PDG values represent an average over 
different data sets and different extraction procedures 
adopted in the various experiments. Consequently, they 
automatically include an enhanced error due to model 
dependence of this extraction. 

The first term in Eq. ( 36 ) is a constant that is model- 



independent, arising from Eq. (34). This model indepen- 
dence reflects the cancelation of the proton-to-resonancc 
transition matrix elements involving the e.m. currents. 
The second term in Eq. (36) , y^, is given by the ratio 



ity amplitudes (we refer the reader to the Appendix for 
details). We summarize the values of ya obtained us- 
ing the PDG values for the helicity amplitudes with the 
respective errors in Table |III[ The lower and upper lim- 
its correspond to taking extreme values of the transition 
helicity amplitudes for the proton and neutron from [40] . 

For the p33(1232) resonance, we assign a conservative 
10% error on its isospin structure. According to the 
PDG, this error should be precisely zero. However, the 
analyzes of Refs. [38, 42] return slightly different results 
for the P33(1232) excitation on the proton and the neu- 
tron, both for real and virtual photons. The discrepancy 
stays below relative 10% for Q 2 < 1 GeV 2 , although this 
conclusion is definitely model-dependent. This observa- 
tion provides motivation for assigning a conservative 10% 
error to yn for the P33(1232). 

Similarly, for the i<37(1950) resonance, the uncertainty 
is driven by the analyses of Refs. [381 H2]- The fit of 
[38] for the proton returns a very mild monopole form 
factor, whereas the neutron data require a dipole form 
factor for the same resonance [32]. Also the strength 
strongly depends on the form of the background, as found 
in our work (see Table [TTJ) . This motivated us to assign a 
conservative 100% uncertainty due to this resonance. 

We note that for both S\\ resonances listed in the Ta- 



ble III the error bar exceeds 100%. This is mostly due 
to the quality of the extracted values for the neutron. It 
is also worth noting that the quark model expectations 
(see Table I of Ref. (TB] for the isospin scaling factors 
within the quark model of Ref. [43] ) are not too far from 
the central values quoted in Table |III| 



1. Uncertainty in isospin rotating the resonances 

To summarize the results of the previous subsection, 
we propose to obtain the contribution of a resonance R 
to the 7Z-interference cross sections cr^ Z £ R by multiply- 
ing the purely electromagnetic cross sections crj 7 £ fl with 
a scaling factor that is independent of W 2 and Q 2 . 
Furthermore, to the precision required here, we rescale 
the trasverse and longitudinal cross sections with the 
same factor. Each such factor contains two parts, as per 
Eq. (36): the first one is model-independent, whereas the 



second one is obtained from the analysis of the proton 
and neutron electromagnetic data, and involves model 
dependence and experimental uncertainties. The values 
of yn are listed in Table |III| with the respective uncer- 
tainties. Correspondingly, for each resonance we simply 
obtain its contribution to the interference structure func- 
tions P7f ,i? from that to the electromagnetic structure 
functions F^ R as 



(37) 



To compute Re OyZ A , we use Eqs. (|20|) and (|16|) with 



of combinations of neutron and proton transverse helic- the input from Eqs. (37) and (25). Finally, we use the 
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parametrizations of the transverse and longitudinal elec- 
tromagnetic cross sections from Model I an d Model II, 
and values of 



sections for photoproduction of vector mesons, 



III 



The uncer- 



factors from Table 
tainty on the contribution of each resonance is obtained 
according to the definition 



AF. 



■yZ,R 
1,2 



(W 2 



Ay R F^ R (W 2 , 



(38) 



where Ay R are the uncertainties quoted in Table |III| Us- 
ing the steps described above for the individual contri- 
butions of resonances to Re O r Za , we can also compute 
the uncertainties A (Re O r Za ) associated with each such 
contribution. Because most resonances do not overlap, 
we treat all these uncertainties as independent, thus we 
define 



A^"ReD 7Z = 



/^|A(ReD« ZA )|a 



(39) 



C. Isospin rotation of the high energy contribution 

We need to employ a well-motivated model to describe 
the isospin dependence of the background contribution. 
One option is to employ the the VDM picture, incorpo- 
rating the simple observation that the photon has the 
same quantum numbers as vector mesons (VM). There- 
fore, it can fluctuate into p, lo or <j> that then scatter off 
the nucleon. This approach underlies the background in 
both Models I and II, so we proceed generally at first. 

According to the VDM, the photon can be represented 
as a superposition of a few vector mesons, 



l7)= E 



h 



■\V), 



(40) 



with fv the VM decay constant. Assuming this basis to 
be complete and orthogonal (no VM mixing), one can 
express the total photoabsorption cross section through 
a combination of total cross sections for vector meson- 
proton scattering, 



E 

V 



fl 



(41) 



At high energies, the total cross section ay P should be 
independent of the VM flavor and the above equation be- 
comes simply a flavor decomposition of the electromag- 
netic total cross section, although this representation is 
of limited use because uy v is unknown. Nevertheless, af- 
ter trivial manipulations this picture leads to the VDM 
(Stodolsky) sum rule |H] that relates the total, real pho- 
toabsorption cross section to a sum of differential cross 



VtotilP) 




16tt — t. [t = 0).(42) 



dt 



This sum rule is based on the assumptions of vector 
meson dominance and almost purely imaginary phase of 
scattering amplitudes at high energy. 

In the naive GVD approach (Model II), this sum rule 
holds only approximately (HERA data: «80%) [45]. The 
missing strength can be attributed to the neglect of non- 
diagonal vector meson-nucleon scattering Vp — > V'p. 
One can then generalize the VDM by including such con- 
tributions by writing down a dispersion relation over the 
vector meson masses. We will denote this non-diagonal 
"continuum" contribution as "X" in the sum over vector 
mesons V. 

Alternately, in the GVD/CDP approach (Model I), in- 
stead of hadronic VM states, the photon hadronic wave 
function (WF) is described in terms of perturbative te- 
states with J = 1. This qq-p&ir forms a color dipole that 
interacts with the target through gluon exchanges. 

Both the naive GVD and GVD/CDP approaches are 
similar in the following instances: they consider the 
interaction of the hadron-like photon with the target 
(hadronic WF for naive GVD, and perturbative qq for 
GVD/CDP), and the interaction of the hadronic states 
is independent of flavor (either VM or quark). This al- 
lows us to cast the ratio of inclusive virtual photon and 
7Z-interference cross sections in the following form: 



(7 



7*P 



E r^ v (W 2 ,Q 2 )a Vp , 



v= P °,uj,4>,x 

a i'Z P = E r i-zv{W 2 ,Q 2 )a Vp , (43) 
v= P °,ui,4>,x 

where "AT" denotes the non-diagonal contribution. 

According to the assumptions of both approaches, 
the flavor factors r 1 * v (W 2 , Q 2 ) and rj* Z v(W 2 ,Q 2 ) only 
contain the information about the projectile (virtual pho- 
ton or Z) and not about the target; this means that they 
cannot depend on the energy but only on Q 2 (the only 
Lorentz scalar that can be constructed from the 7 four- 
momentum) and the flavor of the VM state. However, if 
these flavor factors indeed depend on energy, this would 
signal the breakdown of the models, and would be a 
source of an additional theory uncertainty. For complete- 
ness, we will keep the W 2 dependence. The interference 
flavor factors r^-zv obtain from the purely electromag- 
netic ones using the conservation of the vector current 
(CVC), 
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=0 




e I= 


=0' 






9v 






r 7 



r rZ4> (W 2 ,Q 2 ) = 9 ^r r4> (W 2 ,Q 2 ) = (3-4sm 2 w )r 1 ,4W 2 ,Q 2 ), 



(44) 



for the light flavors. 

With these definitions, we obtain our master formula for rescaling the background contribution: 
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7* p^rZp 

T,L 
7*p->7*p 



(2 - 4sin 2 e w )r^ L Zp (W 2 ,Q 2 ) - 4sin 2 6 w r^(W 2 , Q 2 ) + (3 - 4sin 2 9 w )r^{W\ Q 2 ) + v^^W 2 , Q 



(2 -4 sin 2 



ip v 

,2 , 



rfii(W 2 ,Q 2 ) + <4(W 2 , Q 2 ) + r T <» L x {W\ Q 2 ) 



w ) - 4sin 2 ^^(W 2 , Q 2 ) + (3 - 4 sin 2 w )i# L (W 2 , Q 2 ) 

p p 



T7Z 



1 + R^{W 2 , Q 2 ) + R^(W 2 , Q 2 ) + R^(W 2 , Q 2 

» p p 

I 



(45) 



The ratios R^ L are defined as ratios of transverse (T) 

p 

or longintudinal (L) vector meson (V) production cross 
sections 



The terms 



R 



T,L 



T.L 



l*P^Vp 
7 T,L 



(46) 



'T.L 



T.L 



r 1 *X' T i*ZX 



account for the possible in- 



completeness of the VDM (or three light flavor) basis. 

T L 

For the naive GVD model, Ry_ are obtained from 

p 

the experimentally measured constants fy of the leptonic 
decay V — > e+e~ . Additionally, the presence of the VM 
propagator leads to a prediction for the Q 2 -dependence 

( ™v V 



of each flavor channel 

~f*p^Vp 



Rv — Rv — 



a > 



(jI'p^pp 



f2 4 

J P m v 
f V m P 




T L 

with V — w,(j>. The remaining piece, R x ' is identified 

p 

with the continuum (V — V mixing) contribution. From 



the comparison of the left and right hand sides of the 
VDM sum rule |45j and suplementing this contribution 
with a simple Q 2 dependence to describe the virtual pho- 
toabsorption data at low and moderate Q 2 , one obtains 
i.e. for the transverse ratio 139 



Ri 



= Ro(0) 



(1 + Q 2 /m 2 ) 2 
(l + Q 2 /m 2 ) 



(48) 



with Rc(0) = and m ~ 1.5 GeV. We note that 

due to the monopole Q 2 -dependence of the continuum 
contribution, rather than dipole for the p , the impact 
of the continuum part increases with growing Q 2 . The 
master formula of Eq. ( |45|) - together with the model 
input of Eqs. ((47} and ([48} (see Ref. 3!). for all the 
details of the model) - defines our prescription for the 
isospin rotation of the background contribution within 
the naive GVD model (Model II). 

For pQCD inspired models, such as the GVD/CDP 
used in Model I, the relative strength of the isospin (fla- 
vor) channels is directly related to the quark electric 
charges and is independent of energy and Q 2 . 



T 7->P 



7->w 



7->0 



7-S- J ftp 



(q 



I=0\2 



r 



(q 



s\2 



(<z c ) 



c\2 



((? /=l)2 ' ((? I=1)2 ■ (,1=1)3 



1 : 1 : 2 
9 9 



(49) 



One possible way is to identify the X state in the master 
formula with the cc state, i.e. J/ip. In that case, the 
X contribution in the numerator of Eq. ( 45 ) is given 



according to the SM 

T,L 
r y*ZX _ 
T,L 
IP 



3-. 



sm 

~2~~ 



7 W r>T,L 



(50) 
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The choice of identifying X with J /tp is justified in HERA 
kinematics but is probably less convincing at lower en- 
ergies and low Q 2 . Moreover, the choice X — J/tp and 
the relative strength of different contributions according 
to Eq. (49), corresponds to the VDM sum rule being 



saturated to only 60%, rather than the measured 80%, 
suggesting that it is not very realistic. 

Either way, for the rescaling of the background con- 
tribution in the GVD/CDP model (Model I), Eqs. @ 
and (49) simply combine to a constant factor. Its value 



when using only the three light flavors amounts to 



_7*p->Zp 



, Model I 



2-4 sin 2 9 



w 



1.05 



(51) 



However, when including the charm contirbution, one 
obtains 



(J 7*p-»-7*p 



Model I 



u,d,s,c 



9 



4sin 2 6W ~ 0.85 (52) 



For comparison, a typical value of this ratio within the 
naive GVD Model II (we quote its value at Q 2 = for 
definiteness: in Model II it is Q 2 -dependent, although 
mildly) is 



Model II 



1.92 -4sin 2 6>w« 0.97, (53) 



and a very similar number when including the continuum 
and assuming its size for the jZ cross section to be equal 
to that for the purely electromagnetic case. However, any 
such estimate bears at least 20% uncertainty due to the 
incompleteness of the naive VDM basis and due to the 
unknown flavor structure of the continuum contribution. 

To illustrate the difference in the Q 2 -dependence of the 
total cross section as calculated in Model I and Model II, 
we define the following two ratios, 



i? 77 (IV 2 ,Q 2 ) 



R lZ (W 2 ,Q 2 ) 



Model I 



7 *p->7*p 7 »p_ >7 » p 
T t a L 



Model II 



7*p-^Zp , 7*p->Zp 
°T "L 



~l*P^Zp , 7*P->Zp 
T t a L 



Model I 



Model II ' 



(54) 

where we suppressed the argu ment s of the cross sec- 
tions for compactness. In Fig. 12 we display the Q 2 - 
dependence of i? 77 and R 7 z at two values of W 2 . The 
ratios show very mild JV 2 -dependence, in accord with 
general assumptions used in VDM and GVD/CDP mod- 
els. The Q 2 -dependence shows slight oscillations (at the 
level of 3%) at Q 2 < 2 GcV 2 ; at higher values of Q 2 both 
ratios decrease monotonically, as a result of the naive 
VDM model (Model II) overshooting high-Q 2 data sig- 
nificantly, while GVD / CDP Model I describes data in a 
wide kinematical range. 



-- W 2 = 50GeV 2 
W 2 = 5 GeV 2 




2 3 
Q 2 (GeV 2 ) 



FIG. 12: Ratios R~ n (W 2 ,Q 2 ) (upper panel) and 
R~,z{W 2 ,Q 2 ) (lower panel) are shown as function of 
Q 2 at W 2 = 5 GeV 2 (solid lines) and at IV 2 = 50 GeV 2 
(dashed lines). See text for further details. 



1. Uncertainty in isospin rotating the background 

We wrote our master formula in terms of ratios of 
meson production cross sections, rather than cross sec- 
tions themselves. These ratios were recently measured at 
HERA. The predictions of Eqs. ( |47|49 ) are confronted 
with the experimental data of Ref! [U] at high ener- 
gies and for Q 2 that ranged from zero to several GeV 2 
in Fig. |13| To estimate uncertainties in isospin rotation 
of Models I and II, we will directly compare the model 
predictions of the isospin ratios to the HERA data. The 
common feature of the two models is that these ratios are 
W 2 -independent. Furthermore, Model I predicts them to 
be Q 2 -independent, too. Instead, Model II (naive GVD) 
predicts the Q 2 -running of these ratios. In both cases, 
we will assume that the uncertainty in isospin scaling 
the I — 1 channel (i.e., the p°) is zero. For each flavor 
channel, we define the uncertainty as the discrepancy 



A- "™ 



0-7-S-P 



(Q ) 



a ~1^P 



Model 



V \ ex P 



0-7-S-P 



(55) 



with V = p, oj for VDM and V = p,ui, (J/ip) for pQCD. 
Additionally, for VDM we assign a 100% uncertainty 
to the continuum contribution whose flavor content is 
not defined in the naive GVD approach. Similarly, for 
GVD / CDP model we assign a conservative 100% uncer- 
tainty to the cc contribution, in view of an unsatisfactory 
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FIG. 13: (Color online) Experimental data for ratios of total 
cross sections for elastic vector meson electroproduction in 
comparison with the naive VDM (solid lines) and perturbative 
SU(4) (dashed lines) predictions. 



description of the data for (J/tjj)/p by the 5i7(4) predic- 
tion. 

Arriving at the estimate of the total uncertainty due to 
the isospin structure of the background requires following 
steps. For Models I and II, we insert the uncertainties 



defined in Eq. ( 55 ) and below for each flavor V into the 
master formula Eq. ( 45 1 . Subsequently, we use the rela- 



tion of Eq. (25) and obtain the total uncertainty of the 
interference structure functions, AF^f . We evaluate the 



imaginary part of the dispersion correction of Eq. ( 16 ) 



with AF^- The final step involves evaluating the dis- 



persion integral thereof, Eq. (20). These steps give us 
the uncertainties due to the isospin structure of the back- 
ground within Model I and Model II. To be conservative, 
we will choose the larger of the two as our estimate of 
the non-resonant model uncertainty. 

Anticipating the discussion in the next Section, we note 
that the overall uncertainty is dominated by the contin- 
uum contribution ("AT") within the naive GVD model. 
The only significant assumption about the continuum 
contribution here is that its size (relative to diagonal vec- 
tor meson contributions) is energy-independent , and we 
take it from the data at very high W 2 . Until now, the 
only dedicated study of the VDM sum rule was performed 
at W > 70 GeV at HERA - far from the kinematic re- 
gion that dominates the dispersion integral for 0~ ( z- It 
is not a priori clear that the decomposition of the virtual 
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FIG. 14: (Color online) Sum of resonance and background 
contributions to Re D 7 z for Models I and II as indicated in 
the legend of the plot. 



photon into the VM basis works any differently for 5 GeV 
photons than for 80 GeV photons. 

As part of a program of future measurements to con- 
strain the uncertainties in the dispersion correction, it 
would be useful to have direct data on this sum rule at 
lower energies: 2 < W < 10 GeV. In case that new 
data on the VDM sum rule at these energies will become 
available, it will then be straightforward to include addi- 
tional TU-dependent form factors in Eq. ( 45 1 . Data on 



the virual vector meson photoproduction cross sections in 
this kinematic regime could also provide additional im- 
portant constraints. Together with direct measurements 
fo the inelastic PV asymmetries at these kinematics, such 
measurements could in principle lead to a significant re- 
duction in the quoted theoretical error bar. 



V. RESULTS FOR RED 7 z 

We are now in the position to present results for Q^^in 
the forward direction using the sum ruleofEqs. (U6|20l, 
the Models I and II for the electromagnetic cross sec- 
tions along with the isospin considerations provided in 
the previous sections. We display the sum of resonance 
and background in Fig. [14] 

In Fig. |15[ we display the contributions of various kine- 
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matic regions to Re O^z- The upper panel of Fig, 
evidences that the resonance contribution is dominatecl 
by values of Q 2 < 1 GeV 2 , whereas for the total cor- 
rection, values of Q 2 up to 3 GeV 2 have to be taken to 
saturate the dispersion correction to « 90%. The lower 
panel of that Figure demonstrates that values of W 2 up 
to 25 GeV 2 have to be included under the integration 
to saturate the dispersion correction Re 7 z- The data 
from the resonance region W 2 < 5 GeV 2 (resonance plus 
background) only contribute about 65% of the total. The 
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^33(1232) 


5n(1535) 


O 13 (1520) 


Sn(1665) 


F 15 (1680) 


Pn(1440) 


F 37 (1950) 


£ Res. 


Model I (xl(T 3 ) 


(1.21 ± 0.12) 


(0.28±°;f 7 ) 


(0.18 ±0.03) 


(0.06l°;J 4 6 ) 


(o.o4ir i?) 


(0.09 ±0.03) 


(0.48 ± 0.44) 


(2.34l° ;f ) 


Model II (xl(T 3 ) 


(1.23 ±0.12) 


(0.29t° 3 *) 


(0.18 ±0.03) 


(0-06l ;^) 


(0.041°°!?) 


(0.06 ±0.02) 


(0.40 ±0.36) 


(2.24± ;I1) 



TABLE V: Resonances contributions to the dispersion correc- 
tion to the weak charge of the proton ReD 7 z at the QWEAK 
energy Ei a b = 1.165 GeV, in units of 10~ 3 . For each con- 
tribution, we indicate the uncertainty discussed in the text. 
Results for Model I and Model II are shown. 
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FIG. 15: (Color online) Contributions of different kinematic 
regions to D 7 z for Models I. In the upper panel, contributions 
from various Q 2 ranges are shown: solid curve (full result), 
dashed line (Q 2 < 3 GeV 2 for resonance + background), dot- 
ted line (Q 2 < 1 GeV 2 for resonance ± background). For 
comparison, in the upper panel we display result of integra- 
tion of the resonance contribution to D 7 z over all values of Q 2 
(dash-dotted curve) and Q 2 < 1 GeV 2 (dash-double-dotted 
curve). In the lower panel, contributions to D 7 ^ are shown 
that come from W 2 < 5 GeV 2 (dash-double-dotted curve), 
W 2 < 10 GeV 2 (dash-dotted curve), W 2 < 16 GeV 2 (dotted 
curve), W 2 < 25 GeV 2 (dashed curve), and full result (solid 
curve). 



only over those 



fjj), respectively. 



double integral for Im D 7 z in Eq. ( 16 1 
values of Q 2 (W 2 ) that lie below QaT}' 
After that, the dispersion integral of Eq. (20) is evalu- 
ated without further modifications. 

In Table |IV| we display the background contribution 
as calculated in Model I and II for the QWEAK kinemat- 
ics. It can be seen that the background represents both 
the largest contribution and the source of the largest un- 
certainty. Most notably, within the naive GVD approach 
(Model II) , it is completely dominated by the continuum 
contribution whose isospin structure is undetermined. In 
the pQCD approach (Model I), a contribution similar in 
strength is assigned to the cc state. However, because in 
this case we know exactly how the weak boson couples 
to c-quarks, the uncertainty is about half the size of that 
for Model II. This 50% reduction is simply due to the 
fact that g v ~ \ = \e c . 

The individual resonance contributions are displayed 
in Table [V] It can be seen that the overall uncertainty 





Background 


Model I 


(2.85 ±0.85) xl0~ 3 


Model II 


(3.49 ± 1.92) xl0~ 3 



TABLE IV: Background contribution to the dispersion correc- 
tion to the weak charge of the proton ReD 7 z at the QWEAK 
energy E uh = 1.165 GeV. Results for Model I and Model II 
are shown. 



certainty in two contributions, namely 5n(1535) and 
i 7 37(1950). The former, in turn is dominated by the 
uncertainty in the neutron transition helicity amplitude. 
The heavy resonance state is not well-determined and 
should be studied in greater detail to decrease the re- 
spective uncertainty for the dispersion correction. 



According to the discussion in the previous Section, 
we plot the result for Re D 1 z and display the error bar 
on this calculation in Fig. [TBJ For the central value, 
we take the average of Model I and Model II, and use 
the difference between this central value and either of 
Model I or II as the uncertainty due to modeling the e.- 



m. data. For the isospin rotation-related uncertainty, 
we calculate the error within each model as discussed 
before, and quote the larger of the two. We summarize 
this section by quoting the result of the forward sum rule 
evaluated within two models as follows: 
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Ren iZ (E = 1.165 GeV,t = 0) = [5.46 ± 0.27 (mod. avg.) ± 1.92 (backgr.)iaso (res.)] x 1CT 3 , (56) 

The first uncertainty is due to averaging over the two models, the second due to uncertainty in isospin rotating the 
background, and the third due to isospin rotation of resonances. A possibility of measuring the proton's weak charge 
at Mainz at a lower energy Ei a b — 180 MeV is under consideration presently [JH]> and we quote our prediction for the 
dispersion 7Z correction and the respective uncertainty for that energy, 

Re U lZ (E = 0.180 GeV,t = 0) = [1.32 ± 0.05 (mod. avg.) ± 0.27 (backgr.)±°;^ (res.)] x 10~ 3 . (57) 



We see that the total uncertainty in Re D 7 z is about six 
times smaller at Ei a b = 180 MeV than at Ei a b — 1.165 
GeV. 



VI. ADDITIONAL t-DEPENDENCE OF 
DISPERSION CORRECTIONS 

In the previous section, we provided an educated esti- 
mate for Re D 7 z in the exact forward direction. How- 
ever, real experiments are carried out at finite momentum 
transfer t, in particular \t\ = 0.03 GeV 2 for the kine- 
matics of the QWEAK experiment. To extrapolate the 
forward sum rule to non-zero momentum transfer, we 
employ the phenomenological model that was success- 
fully used for the beam normal spin (Mott) asymmetry in 
elastic ep-scattering [4"9"H5"T] . This model is inspired by i) 
experimental data on the Compton differential cross sec- 
tion at small t and high energy, and ii) the assumption 
of the predominantly imaginary phase of the Compton 
amplitude at high energies (as for the pomeron). 

The data exhibit an exponential i-dependence, 

dt \dtj t=0 
with the slope parameter B = 7± 1 GeV~ 2 52] . The 



differential cross section is related to the Compton am- 
plitude squared, whereas the total cross section - through 
the optical theorem - is related to the imaginary part of 
the Compton amplitude. Naively, then, one might expect 
the t-dependence of the total cross section near the for- 
ward scattering limit to be close to half as rapid as that 
of the differential cross section. Based on this ansatz, 
Rcf. [49 proposed parameterizing the i-dependence of 
the slightly off-forward total cross section as 

otot(t) « a tot (* = 0) e"^ . (59) 

This parameterization becomes precise at very high ener- 
gies where the cross section is pomeron-dominated. We 
will follow a similar parameterization here. 

This intrinsic i-dependence of the 7.Z-box contribution 
should be combined with the 77-box contribution that 
becomes non-zero when going to finite t. We found the 
effect of the dispersive contributions to d~~ on A pv to 
be negligibly small, of order below 0.1% at the QWEAK 
kinematics at -t = 0.03 GeV 2 . The reason for this small- 
ness is due to an explicit i-suppression of D 77 with re- 
spect to the tree level PC amplitude. Using the same 
approach, we obtain for the i-dependence of the disper- 
sion correction 

D yZ (E,t) = n, z ( El 0) eXp{ -^ /2 \ (60) 
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iQWEAK (t = -0.03 GeV 2 ) 
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-t (GeV 2 ) 

FIG. 17: (Color online) The dependence of the combined dis- 
persion correction to on the elastic momentum transfer t. 
The central value is shown in red; the thin black lines indicate 
the uncertainty. 



according to the definition of \3 1 z as the ratio of the 
7Z-box contribution to the PV amplitude f^E. t) to the 
elastic proton electromagnetic form factor F^ p (t). 



In Fig. [T7J we display the t-dependence of the com- 
bined dispersion correction for small values of the elastic 
momentum transfer. It can be seen that one can expect 
that at \t\ — 0.03 GeV 2 , the dispersion correction de- 
creases by only about 2% relative to its value at t = 0, 
and the same is valid for the uncertainty in calculating 
this correction. We emphasize, however, that the model 
for the t-dependence is derived from high-energy Comp- 
ton data, and it is not necessarily applicable to the reso- 
nance contributions. Thus, our estimate of the effect of 
the t-dependence should be considered as an exploratory 
investigation. Having this caveat in mind and taking 
into account this t-dependence, we obtain our final re- 
sult for the dispersion -fZ correction at the kinematics of 
the QWEAK experiment: 



ReD^z A (E = 1.165 GeV, t = -0.03 GeV 2 



(61) 



5.39 ± 0.27 (mod. avg.) ± 1.88 (backgr.)j;°;g (res.) ± 0.07 (t - dep.) 



x 10" 



To assess the relative impact of the energy-dependent 
contribution from 7 z A we first quote the result from 
Refs. [31 E] for the weak charge, as defined in Eq. 

Q p w = 0.0713 ± 0.0008 (62) 

Compared to this prediction, the relative effect of the 
7 z A contribution at the kinematics of the QWEAK ex- 
periment is 

Re ^ ZA = (7.6 ±2.8)% (63) 
Qw 

Because this contribution was initially neglected in the 
analysis of radiative corrections for the QWEAK exper- 
iment, the final theory prediction and the respective un- 
certainty have to be corrected to include it. Treating 
all the individual uncertainties quoted above as indepen- 
dent, we obtain 

= (0.0767 ± 0.0008 ± 0.0020 7 z). (64) 

As discussed in Refs. [3l [6], the first error includes a 
conservative ~ 1% error associated with C^z that ap- 
pears in the prediction for Q^- The additional uncer- 
tainty associated with the energy-dependent contribution 



from 0-yz A i s almost entirely due to the uncertainty due 
to the isospin structure of the background. The latter, 
in turn, is largely dominated (70%) by the uncertainty 
due to the "continuum" contribution that does not have 
well-defined isospin content. 

To recollect, the continuum contribution that arises in 
both the naive GVD and GVD/CDP frameworks is a 
measure of the incompleteness of the vector meson ba- 
sis for the energetic photon. Its value at the real pho- 
ton point is obtained from the direct comparison of the 
VDM sum rule with the experimental photoproduction 
cross sections where one finds roughly a 20% deficit in the 
naive GVD approach. Departing from the real photon 
point, one employs phenomenological models for the (in- 
dependence for this contribution, such that in this way 
the generalized VDM description fits the virtual photoab- 
sorption data at moderate Q 2 . It turns out that above 
Q 2 = 2 GeV 2 , the continuum contribution becomes dom- 
inant. 

Strictly speaking, these observations only apply at high 
energies, as the VDM sum rule measurement at HERA 
was performed at W = 82 GeV. In absence of an in- 
dependent evaluation at lower eneriges, we are forced 
to extrapolate this isospin decomposition down to lower 
energies. While this extrapolation is in line with the 
general assumptions of the VDM, there is no guarantee 
that the isospin decomposition of the photon wave func- 
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tion is energy- independent. To illustrate where the high- 
energy assumptions may break down, we note that one 
of the purely high energy scattering assumptions in the 
derivation of the VDM sum rule is the neglect of the real 
part of the forward Compton amplitude with respect to 
the imaginary part. While this holds for the Pomeron 
- whose phase is almost purely imaginary - Reggeon 
exchanges contribute to both real and imaginary part. 
The main contribution to ReD 7 ^ comes from energies of 
W < 5 GeV where the reggeon contribution dominates. 
Thus, a re-evaluation of the VDM sum rule at JLab en- 
ergies will likely help to reduce the theory uncertainty on 
the ReD 7 2 calculation. 



VII. COMPARISON TO RELATED WORKS ON 
DISPERSION jZ CORRECTION 

In this Section, we briefly outline the main improve- 
ments achieved in this article with respect to our pre- 
vious work, as well as the recent work carried out by 
other groups. In Ref. |llj . the forward dispersion rela- 
tion for \3 7 z was derived and evaluated with the result 
of ReD 7 z w 0.003. However, that study used an over- 
simplified model of virtual photoabsorption and for the 
sake of simplicity assumed that F^f = F±2- These as- 
sumptions did not allow for a realistic study of uncer- 
tainty of that result. 

In Ref. |16j . we improved on these two points: we em- 
ployed a phenomenological model of Bosted and Christy 
that fits virtual photoabsorption data over a large kine- 
matic range; considered the isospin structure of each con- 
tribution; and discussed the possible ways to estimate un- 
certainty on this calculation. Although in [16] we were 
able to develop the general method that we use in this 
manuscript, no robust theory error bar was obtained. 

Sibirtsev, Melnitchouk, Blunden and Thomas in Ref. 
[T5] re-checked the findings of Ref. [TT] . That group care- 
fully re-derived the sum rule pointing out two errors in 
[llj which have been corrected here. Furthermore, the 
authors of [TS] proposed a model of the virtual photoab- 
sorption that was directly fit to the experimental data. 
This allowed them to obtain an estimate for the uncer- 
tainty in the dispersive calculation of D 7 z from the error 
bar of the fit. To obtain the interference structure func- 
tions from the purely electromagnetic ones F^, the 
authors of |15j relied on isovector dominance in the res- 
onance region. For the background, employed a simple 
scaling prescription for the background 
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extrapolating the isospin structure from the DIS region 
to low energies. Ref. [15] confirmed that the dispersion 
correction is sizeable, obtaining the value quoted in Table 
[TJ where the error bar is due solely to fitting the electro- 
magnetic data. The uncertainty on the isospin rotation of 



the electromagnetic data was not included. Correspond- 
ingly, the error bar quoted in [15] only contains one part 
of the total theoretical uncertainty in evaluating ReO jZ . 

Rislow and Carlson subsequently performed another 
computation of the dispersive contribution to □ 7 z A [IZ] • 
These authors again confirmed the derivation of the sum 
rule and re-evaluated it obtaining the somewhat larger 
result quoted in Table [TJ This was achieved by using a 
different model for the resonances (the same as used in 
our present work) that allows for a better fit of electro- 
magnetic data. The central value of Ke\3 1 z is very close 
to our estimate in forward direction, see Eq. (56 1. The 



background was taken in a phenomenological form and 
continued into the DIS region. The authors discussed in 
some detail the procedure of isospin-rotating the reso- 
nance contributions using the constituent quark model. 
For the background, Ref. [T7] follows to a large extent 
the isovector dominance picture with a perturbation on 
top of that that provides an estimate of the uncertainty. 

Both Refs. |T5l [17] argue that the uncertainty on the 
dispersive calculation of ReD 7 z is well under control and 
can be easily accomodated within the error budget of 
the QWEAK. However, because they do not provide a 
model-independent analysis of the isospin structure, we 
believe that the estimates of the error bars obtained in 
those two works is unlikely to be complete. 

In this manuscript, we believe we have developed the 
most robust and model-independent estimate to date of 
the absolute size and the uncertainty of the dispersion 
*yZ correction to Apy in the forward limit. We used the 
most recent fit to resonance data, supplemented by two 
different models of the background. We demonstrated 
that the two models used in this work indeed provide a 
good description of the experimental data in a very wide 
kinematic range of two variables W, Q 2 . 

For the 1 = 1/2 resonances, we employed an isospin 
rotation that is reminiscent of that for the elastic elec- 
troweak form factors. This allows one to unambiguously 
relate the ratio of interference jZ and the electromag- 
netic cross sections to combinations of transition helic- 
ity amplitudes for the photoexcitation of a given res- 
onance on the proton and the neutron. We used the 
most recent values and uncertainties for the latter from 
PDG 40 . The main sources of the uncertainty for the 
1 = 1/2 resonances is the neutron transition helicity am- 
plitude of Sii(1535), where a more precise extraction of 
the transition helicity amplitude on the neutron would 
be needed. For the 7 = 3/2 resonances, the isospin ro- 
tation is straightforward. However, we assigned a 100% 
uncertainty to the contribution of the heavy resonance 
state that lies close to F 3 7(1950) but cannot be reliably 
identified with the latter and whose isospin structure is, 
therefore, uncertain. 

For the background we utilized two models based on 
the framework of vector dominance model (VDM) that 
provides a prescription for the isospin decomposition of 
total photoabsorption at high energies. The VDM sum 
rule states an equality between the total photoabsorption 



22 



and differential cross sections for forward vector meson 
production. This sum rule has been tested experimen- 
tally, albeit only at very high energies. The precision to 
which this sum rule holds provides us with one handle for 
assessing for the robustness of our isospin decomposition 
of the electromagnetic data. To investigate the model 
dependence, we use two different models that obey the 
general requirements of the VDM but originate from two 
kinematically distinct regimes: Model I is a pQCD color 
dipole model (what we have called the GVD/CDP ap- 
proach) that is continued down to the real photon point 
by employing phenomenological input |37j . Model II is 
a "naive" generalization of VDM quoted in [32. • The 
two models lead to similar numerical results, but within 
each model the estimate of uncertainties is different. The 
largest contribution to our quoted theoretical uncertainty 
arises from lack of knowledge of the isospin structure of 
the terms in these models that are not uniquely associ- 
ated with any one of the three lightest vector mesons. To 
be conservative, we have chosen the largest of the corre- 
sponding uncertainties from the two models. 

Finally, we considered a phenomenological model for 
the intrinsic ^-dependence to extrapolate the forward sum 
rule to the experimental kinematics. We find that the ef- 
fect of such extrapolation is not significant. However, we 
consider this approach to be exploratory, and an addi- 
tional uncertainty on i-dependence may have to be taken 
into account. 

To reduce the uncertainty associated with D 7 z to a 
level below 2%, there exist a number of avenues that 
could be pursued. The most direct would be to per- 
form measurements of the inelastic PV asymmetries in 
the kinematic region that dominates the dispersion inte- 
gral: W < 5 GeV and Q 2 < 3 GeV 2 . Doing so would 
provide information on the electroweak structure func- 
tions that enter the dispersion integral for 
thereby mitigating the need for a model with which to 
carry out the isospin rotation. Additional constraints 
could be obtained by experimentally testing the VDM 
sum rule at the lower energies relevant to the aforemen- 
tioned kinematics; by performing precise measurements 
of the electromagnetic neutron-to-resonance transition 
cross sections, thereby yielding the corresponding helicity 
amplitudes — particularly for the 511(1535); and identi- 
fying the isospin of the ^37(1950) resonance. 

An alternate strategy would be to perform a mea- 
surement of Apv at lower energy, given that the mag- 
nitude of, and uncertainty in, \3 7 z decrease monotoni- 
cally with decreasing energy as indicated in Fig. [16] and 
Eq. (57). From the standpoint of probing physics be- 



yond the Standard Model, a measurements of atomic PV 
observables for different isotopes may also be interest- 
ing. The largest atomic theory uncertainties cancel from 
ratios of these observables [55]. and the leading sensitiv- 
ity to new physics is dominated by the effects on the 
proton weak charge [4]. To the extent that uncertain- 
ties in the neutron distributions can be constrained (e.g., 
through measurements of the elastic PV asymmetry for 



heavy nuclei) , "isotope ratio" experiments may provide a 
cross-check on any inferences about new physics derived 
from the QWEAK measurement. Given the experimen- 
tal and theoretical challenges involved in each of these 
efforts, an ideal program may entail a combination of the 
aforementioned measurements. 
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Appendix A: Isospin rotation of the resonance 
contr ibut ions 



In Standard Model, the Z and 7 couplings to the 
quarks are related by an isospin rotation, 



J em — 1 J 1=0 ^ 1 J 1=1 + y J s 
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The e.m. charges given by q — — ^= 
q s = — |, and the weak charges are g(T° = 



1 , q I=1 
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723 S b W, 



l+is 2 9w Consequently, 



giT 1 = ^(2-4s 2 ^), g» v - , , 
this isospin decomposition relates weak proton form fac- 
tors to the proton and neutron electromagnetic form fac- 
tors, 

Wnc,v\p) = (1 - 4sm 2 8w)WL\p) ~ H-O) 

(A3) 

Above, we neglected strangeness contributions that are 
generally small. 

A similar relation is valid for I = \ resonances, as well: 



(X\J»\p) = (1 - 4sin 2 9w)(X\J*Jp) - (X\J? rr 



(A4) 



Then, the contribution of a resonance R with isospin 
1/2 to the interference jZ cross section can be related to 
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its contribution to the electromagnetic cross section by To proceed, we use the definition of the transition he- 

licity amplitudes, 

<p\JL\R)W$rc,v\P) = (l-is 2 w )\(R\J^\p)\ 2 

- {p\J? m \R)(R\J? m \n) (A5) 



ifi/afl^.Q 2 ) = ( R > A * = 1 /2|J e Al m (A 7 = l)\ P (n),A N = -1/2) 
A P $ /2 (W 2 ,Q 2 ) = (R,A R = 3/2|J£ m (A 7 = l)\ P (n),A N = 1/2) 
S P R %{W 2 ,Q 2 ) = (R,A R = 1/2|J C ^(A 7 = 0)\p(n),A N = 1/2) 



(A6) 



where we introduced photon helicity A 7 = 0, ±1, nucleon 
helicity A^r = ±1/2 and the helicity of the resonance 
R that is related to the former two as Ar = A^r + A 7 . 
Resonance contributions to the total cross sections <jt,l 
arc related to the helicity amplitudes as 
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with Mr,Tr and qr the resonance mass, width and the 
three-momentum of the virtual photon on the resonance 
position, respectively. In the above equation, the argu- 
ments W 2 , Q 2 of the cross sections and helicity ampli- 
tudes were suppressed. 



We combine the definition of Eq. ( 35 ) with Eqs. ( A5 1 , 



(A6), (A7), and finally obtain 
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For spin-1 resonances, only A p 'p 2 pieces contribute in 
the transverse ratios C^/ 7 - To a good approximation, the 
width and position of a resonance can be assumed to be 
the same for proton and neutron induced reactions. In 
this case, the VF-dependence cancels out in the ratio, and 
it is fucntion of Q 2 only. 

We write in general 



e zh {Q 2 ) 

Cf /7 (Q 2 ) 



(1 - As 2 6 w ) - y R x x R (Q 2 ), 
{l-As 2 e w )-y R x x R (Q 2 ), 



(A9) 



with yR ({jr) the values of the ratio of the neutron and 
proto n tr ansverse (longitudinal) helicity amplitudes in 
Eq. (A6) at Q 2 = 0, and xr (xr) the respective form 



factors. The form factors are normalized to unity at the 
real photon point. 



For the resonances of isospin 3/2, the transition is 
purely isovector, and the ratio of the cross sections is 

given by = 2 — Asin 2 9 w , and is Q 2 -independent. 
However, for the phenomenological analyses of the in- 
clusive virtual photoabsorption data on the proton and 
neutron, Refs. [351 S2] widely used in this work, this 
rule does not hold. For the A(1232) it holds to about 
10%. For the F 3 7(1950), the proton and neutron transi- 
tion form factors show very different behavior (monopole 
for the proton vs. dipole for the neutron). Furthermore, 
the unnaturally mild monopole form factor raises a ques- 
tion of whether this contribution should be considered 
as part of the background where monopole form factors 
arise naturally in the VDM picture. Correspondingly, 
rather than operate with a form factor xr(Q 2 for the 
two isospin-3/2 resonances we will assign an uncertainty 
to the ratios y R =3/2 = -1: 10% for the A(1232) and 
100% for the F 37 (1950), and use x R {Q 2 ) = 1 for both. 

We next turn to the form factors xr(Q 2 ) of the 
isospin-1/2 resonances. To estimate these, one needs 
the Q 2 -dependence of the transition helicity amplitudes 
for the excitation of these resonances. Unfortunately, 
the phenomenological fits of Refs. [3H H2] do not 
provide us with this information: they only give us 



A p A' 



|A^™y 2 | 2 and |-5^"L| 2 . Instead, we need e.g., 
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For spin-1/2 resonances, only 

A 1 / 2 's contribute. Then, one has for the second terms in 
Eq. ((A8j) 
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and the only missing piece above is the relative sign of 
the proton and neutron helicity amplitudes. This sign 
is well defined and can be taken, for instance, from the 
PDG or from quark model [43]. For spin-3/2 resonances 
Di 3 (1520) and i*\ 5 (1680), the information provided by 
Refs. [351 S2] is not sufficient to determine respective 
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xr(Q 2 )- We can only approximately estimate those by 
noticing that for these resonances, the PDG suggests that 
the p — > N* transition is completely dominated by the 
A 3 / 2 helicity amplitude 40 . Then, we can adapt the 
same logic as for the spin-1/2 resonances, by substituting 



A 3 / 2 y s in place of Ai/ 2 's in Eq. (A10) 



However, this procedure cannot be considered reliable 
since such "extracted" form factors xr(Q 2 ) will contain a 
model dependence that is very hard to estimate. Instead, 
we will use the following reasoning. We verified that with 
the approximations described above, the results of Refs. 
[551 |4"2"] lead to the form factors xr(Q 2 ) that differ from 
1 by at most 10-20% for values of Q 2 < 0.6 - 0.8 GeV 2 
for all five isospin-1/2 resonances. At the same time, the 
PDG quotes the errors for the helicity amplitudes [30] for 
the excitation of those resonances that are conservative 
enough to accomodate these 10-20% discrepancy. Indeed, 
the PDG values represent an average over world data and 
over various analyses, therefore the errors that they quote 
contain not only the statistical and systematic error of 
each experiment, but also the systematic error due to 
model dependence of those analyses. This means that at 
low values of Q 2 , the error introduced if setting xr(Q 2 ) ~ 
£r(0) = 1 is reasonably small as compared to the error 
in the respective i/r. While at larger values of Q 2 this is 
not the case any longer, due to resonance form factors the 
impact of these values of Q 2 on the dispersion correction 
□ 7 z is small. From the discussion of the results, we see 
that 

• the resonance contribution is dominated by the 
A(1232) for which the issue of the uncertainty in 



xr(Q 2 ) is controlled within 10%, as discussed ear- 
lier; 

• the overall uncertainty on the resonance contribu- 
tion is dominated by that due to the problem of 
the identification of the high lying resonance in the 
analysis of Bosted and Christy with the ^37(1950); 

• the total uncertainty in the dispersion correction 
□ 7 Z is dominated by the uncertainty due to the 
background contribution. Then, even doubling the 
uncertainty in the contribution of the Sii(1535) 
due to xr(Q 2 ) will not significatly change our over- 
all conclusions. 

This allows us to set all xr(Q 2 ) = 1 for all seven res- 
onances considered here (including the isospin-3/2 reso- 
nances discussed earlier). The error introduced by this 
approximation is safely covered by using the conservative 
PDG errors for the resonance helicity amplitudes. 

Finally, we discuss the ratios of the longitudinal cross 
sections Cz/ 7 - I n R- e f- [32]) it was shown that the hypoth- 
esis that the ratio of the resonance contributions to lon- 
gitudinal and the transverse cross sections for the proton 
and for the neutron target are equal, is well supported 
by the experimental data. Although this conclusion is 
model-dependent, as well, the general impact of the lon- 
gitudinal cross section on the \3 7 z was found to be very 
small. This allows us to use the assumption of Ref. [42 
here and set Cz/j = £z/i f° r an seven resonances. As a 



result, we arrive at Eq. (|37|) with 



1-4 sin" 9i 



,zh — *~ "* Din Jw-Vr 
where the values of y# are listed in Table |III| 
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